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ABSTRACT 

We develop a general expression for the 7 — 7 absorption coefficient, a 77 , for 7-rays propa- 
gating in an arbitrary direction at an arbitrary point in space above an X-ray emitting accretion 
disk. The X-ray intensity is assumed to vary as a power law in energy and radius between the 
outer disk radius, Ro, and the inner radius, R ms , which is the radius of marginal stability for a 
Schwarzschild black hole. We use our result for a 77 to calculate the 7 — 7 optical depth, r 77 , for 
7-rays created at height z and propagating at angle $ relative to the disk axis, and we show that 
for $ = and z > Rq, t 77 oc E a z~ 2a ~ 3 , where a is the X-ray spectral index and E is the 7-ray 
energy. As an application, we use our formalism to compute the minimum distance between the 
central black hole and the site of production of the 7-rays detected by EGRET during the June 
1991 flare of 3C 279. In order to obtain an upper limit, we assume that all of the X-rays observed 
contemporaneously by Ginga were emitted by the disk. Our results suggest that the observed 
7-rays may have originated within < 45 GM/c 2 from a black hole of mass > 10 9 M Q , perhaps in 
active plasma located above the central funnel of the accretion disk. This raises the possibility of 
establishing a direct connection between the production of the observed 7-rays and the accretion 
of material onto the black hole. We also consider the variation of the optical depth as a function of 
the angle of propagation <I>. Our results indicate that the "focusing" of the 7-rays along the disk 
axis due to pair production is strong enough to explain the observed degree of alignment in blazar 
sources. If the 7-rays are produced isotropically in 7-ray blazars, then these objects should appear 
as bright MeV sources when viewed along off- axis lines of sight. 
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1. INTRODUCTION 



Over the past several years, the EGRET instrument on board the Compton Gamma Ray 
Observatory (CGRO) has detected high-energy 7-rays from 44 active galactic nuclei (AGNs), and 
several of these objects have also been detected at lower energies by COMPTEL and OSSE (Fichtel 
et al. 1994; Hartman et al. 1992; Hermsen et al. 1993; McNaron-Brown et al. 1994; von Montigny 
et al. 1995). One of the most dramatic events seen was the intense 7-ray flare of 3C 279 observed in 
June 1991, which has called into question the viability of a number of theoretical models for 7-ray 
emission in AGNs. Most of the energy emitted during the flare appeared at energies exceeding ~ 
100 MeV, with an implied (isotropic) luminosity of ~ 10 48 ergss -1 . However, the actual luminosity 
may have been several orders of magnitude lower if the emission was not isotropic, but rather 
"focused" onto a small region of the sky due to some alignment mechanism such as relativistic 
beaming or geometrical collimation. In particular, the hypothesis of relativistic motion is consistent 
with observations of apparent superluminal motion, flat radio spectra, rapid optical variability, and 
high polarization in blazars. 

In addition to beaming and geometrical collimation, 7 — 7 pair production in the X-ray field 
of the accretion disk could in principle lead to preferential escape of the 7-rays along the symmetry 
axis of the disk, due to the strong angular dependence of the pair-production cross section. The 
phenomenon of 7 — 7 "focusing" is related to the more general issue of 7 — 7 transparency, which 
sets a minimum distance between the central black hole and the site of 7-ray production, z 77 . One 
usually finds that z 11 greatly exceeds the cooling length, 



for relativistic electrons of energy E traveling at essentially the speed of light c and suffering inverse- 
Compton and synchrotron losses in the intense magnetic and radiation fields close to the center of 
the active nucleus (Dermer & Schlickeiser 1993). An efficient pair creation and/or reacceleration 
mechanism is therefore required in order to explain the presence of relativistic electrons sufficiently 
far from the black hole to produce observable 7-rays. The nature of this mechanism is unclear, and 
a variety of possibilities have been suggested, including shocks propagating along a jet (Blandford 




(1.1) 
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1993), the decay of charged pions resulting from hadronic collisions (Bednarek 1993), and the ac- 
celeration of electrons and positrons via wave-particle interactions in the funnel of a thick accretion 
disk (Becker, Kafatos, & Maisack 1994). Constraints on the mechanism are provided by estimates 
of the distance over which it must operate, z 77 — -Zcooi' 

Clearly, the closer one can place the 7-ray 
production site to the black hole while maintaining 7-ray transparency, the easier it will be to 
power the observed emission using accretion. It is therefore interesting to calculate the minimum 
distance satisfying the 7-ray transparency constraint in the context of a specific X-ray emission 
scenario, and that is the focus of this paper. 

Gould & Schreder (1967) calculated the absorption coefficient (a 77 oc cm -1 ) describing the 
process 7 + 7' — > e + + e~ for high-energy photons traversing a soft photon gas, and their formalism 
has been used subsequently by Blandford (1993), Dermer & Schlickeiser (1994), and Bednarek 
(1993) to estimate the absorption probability for 7-rays traversing the X-ray field of an AGN. The 
detailed results depend strongly on the energy and angular distribution of the soft radiation, and are 
therefore sensitive to the presence of scattered radiation, which is more isotropic than the primary 
radiation emitted by the disk. Blandford (1993) and Dermer Sz Schlickeiser (1994) assumed that the 
soft photon distribution is isotropic, which is reasonable if the X-ray distribution is dominated by 
the scattered component. However, the relative importance of the scattered photons compared to 
the primary (accretion disk) photons is determined by the electron scattering optical depth above 
the disk, which is not known with any certainty. Bednarek (1993) adopted the opposite viewpoint 
in his calculation of the optical depth to pair production (r 77 ), neglecting electron scattering and 
instead treating the full anisotropy of the soft radiation field produced by the accretion disk. He 
assumed a standard disk radiating as a blackbody with local temperature T(R) oc i?~ 3 / 4 . The 
resulting X-ray intensity (integrated over the face of the disk) has the characteristic power-law 
shape I v oc z/ 1 / 3 (e.g., Pringle 1981). This is inconsistent with most AGN (and blazar) X-ray 
spectra (with I v oc v~ x l 2 typically), which are usually better fit by hot, two-temperature accretion 
disk spectra (e.g., Shapiro, Lightman, & Eardley 1976, hereafter SLE) in the energy range relevant 
for pair production (e > 261 eV for 7-ray energy E = 1 GeV; see eq. [3.4]). 

In 7-ray loud blazars such as 3C 279, variability arguments suggest that a significant fraction 
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of the observed X-rays are probably produced in a relativistic beam moving toward the observer 
(Makino et al. 1989). Since these X-rays are coaligned with the observed high-energy 7-rays, they 
do not contribute substantially to the 7 — 7 opacity, and therefore they do not strongly constrain 
the minimum distance between the black hole and the site of 7-ray production. However, as we 
demonstrate below, an estimate of this distance can be obtained by considering the attenuation of 
the 7-rays due to collisions with X-rays emitted by an underlying accretion disk. In order to tie the 
distance estimate we obtain to the multifrequency observations of particular objects such as 3C 279 
without incurring additional theoretical complexity, in this paper we make the further simplifying 
assumption that all of the contemporaneously observed X-rays are emitted by the accretion disk. 
By overestimating the X-ray flux emitted by the disk, we obtain an "upper limit" on the minimum 
distance between the black hole and the point of creation of the observed 7-rays, and we expect 
that a proper accounting for the component of the X-ray flux generated in the relativistic beam 
would further reduce this distance. In any event, the upper limit we obtain is small enough to 
allow for the possibility of establishing a direct connection between the production of the observed 
high-energy 7-rays and the accretion of matter onto the black hole. 



Electrons in the hot, inner region of a two-temperature accretion disk produce power- 
law inverse-Compton X-ray emission by upscattering UV radiation produced in the cool (single- 
temperature) surrounding region. The temperature of the electrons in the hot region is maintained 
at a nearly uniform value in the range T e ~ 10 9 — 10 10 K by the "Compton thermostat" mechanism 
(SLE; Eilek & Kafatos 1983, hereafter EK), and the emitted X-ray spectrum has the characteristic 
power-law shape 



in the frame of the source, where y = (4 k T e /m e c 2 ) Max (r es , Tg S ) ~ 1 is the Compton y-parameter, 
m e is the electron mass, and r es is the electron-scattering half-thickness of the disk. The power-law 
extends up to photon energies e ~ kT e , above which the spectrum declines exponentially. Since a 
and T e are essentially constant in the hot inner region of the disk, it follows that the radial variation 



2. X-RAY INTENSITY OF THE DISK 
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(2.1) 



6 



of the X-ray intensity follows that of the total flux, given by 



3 GMM 
8^ R 3 



6GM 
~&R 



) 



1/21 



F(R) = 



1 - 



(2.2) 



for a quasi-Keplerian accretion disk around a Schwarzschild black hole of mass M, where M is 



the accretion rate (Shakura & Sunyaev 1973). Note that the total flux has a nearly power-law 
dependence on radius F(R) oc R~ 3 , with a sharp cutoff at the radius of marginal stability for 
circular orbits, R ms = 6R g , where R g = GM/c 2 . 

Based upon equations (2.1) and (2.2), we conclude that the X-ray intensity in the frame of 
the disk can be represented approximately by 



where a is given by equation (2.1), R is the outer radius of the two-temperature region, and 
Rms < R < Rq- According to equation (2.2), lo = 3 for a two-temperature disk, but we shall leave 
w as a free parameter so that we can make comparisons with uniformly-bright disks, with to = 0. 
Since it is unclear to what degree electron scattering isotropizes the X-rays above the disk, we follow 
Bednarek's approach, and neglect the effect of electron scattering on the angular distribution of the 
X-rays. Our neglect of a scattered X-ray component is in some sense self-consistent, since it turns 
out that most of the 7 — 7 attenuation occurs very close to the disk, where the X-ray distribution 
is thought to be dominated by unscattered accretion disk photons (Dermer &: Schlickeiser 1994). 
However, the scattered radiation may attenuate the 7-rays over much larger length scales. 

Although the inner region of the disk terminates at the radius of marginal stability, the 
outer radius R$ of the two-temperature region is more difficult to quantify because of its depen- 
dence on the value (and possible variation) of the viscosity parameter a, but one typically finds 
30i? 9 < Rq < 100 R g (SLE; EK). Most of the photons emitted in the "cool" region beyond radius 
i?o are too soft to produce pairs via collisions with 7-rays of energy E ~ 1 GeV, and therefore the 
effect of these photons on the 7 — 7 absorption coefficient is negligible. 

In order to link our calculation of the 7 — 7 absorption coefficient with the X-ray and 7- 
ray data, we must express the disk-frame intensity parameter Iq in equation (2.3) in terms of the 
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parameter Fq characterizing the observed (Doppler-shifted) X-ray flux, 
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m e c z 



which is related to the observed intensity 



I' e {z',R)=I' 



(oc s cm ) 



R_ 

Rg 



(2.4) 



(2.5) 



by 



F[ = j I' e cos9' dW , (2.6) 

where primes denote quantities measured in the frame of the earth. 

Since we are mainly interested in blazar sources, where we are presumably viewing the 
accretion disks face-on, we can express Fq in terms of I' Q and the luminosity distance to the source 
D by substituting for the intensity in equation (2.6) using equation (2.5) and integrating over the 
face of the disk. Making the substitutions dQ,' = 2ir sin 0' d9' and tan#' = R/D yields 



Rq = 2^D 2 / / -fr 
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or, in the limit D R , 
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for u ^ 2, where I' Q is related to Iq via the Lorentz transformation of the intensity, which gives for 
a source located at redshift Z (e.g., Begelman, Blandford, & Rees 1984) 



I' = {l + Z) 



-3 — a j 



(2.9) 



In § 3 we obtain a general formula for the 7 — 7 absorption coefficient for 7-rays propagating above 
an accretion disk radiating with the X-ray intensity distribution given by equation (2.3). 

3. PHOTON-PHOTON PAIR PRODUCTION 



The exact cross section for 7 — 7 pair production is given by the well-known result (e.g., 
Jauch & Rohrlich 1955), 



a 11 {l3) = -a T {l-(3 2 ) 



(3-/3 4 )ln(i±|)-2/?(2-/3 2 ) 



(3.1) 



where (3 is the velocity of the positron (or electron) in the center-of-momentum (CM) frame in units 
of c, and a T is the Thomson cross section. The Lorentz factor of the particles in the CM frame, T, 
is related to the energy of the soft photon e, the energy of the 7-ray E, and the interaction angle 
between the photon trajectories 9 by 

r2 _ 1 _ e£(l-cosfl) 

where e, E, and 9 are all measured in the same (arbitrary) frame. The cross section for the process, 
cr 77 , vanishes as T — > 1. Hence for fixed e and E, the absorption probability decreases as 9 — > 0, 
and consequently the optical depth to pair production (r 77 ) drops sharply as the photons become 
co-aligned. We therefore anticipate that the value of r 77 will depend sensitively on the angular 
distribution of the soft radiation relative to the direction of propagation of the 7-ray. 

3.1. 7 — 7 Absorption Coefficient 

Working in terms of the specific intensity I e , we can write the differential contribution to 
the absorption coefficient a 77 due to soft photons propagating in energy range de and in directional 
range dQ as 

da 11 = — cr 77 ((3) (1 — cos 9) dedtt . (3.3) 

We assume that the X-ray intensity of the disk is given by equation (2.3), and that any change in 
the X-ray spectral index a occurs below the minimum soft-photon energy required to produce a 
pair with T = 1, 

Taking E = 1 GeV, we find that e m ; n = 0.522 keV for a right-angle collision (cos 9 = 0), and 
e m in = 0.261 keV for a head-on collision (cos 9 = —1), which is the minimum possible X-ray energy 
resulting in pair production. The X-ray spectrum given by equation (2.3) would conflict with the 
measured EGRET spectrum if extended to arbitrarily high energies, and consequently a break must 
exist in the energy range 0.1 — 1 GeV. However, for our purposes, the X-ray power-law can formally 
be extended to infinite energy without producing numerical divergence because cr 77 oc e _1 lne as 
e — > 00 for fixed E and 9. 
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Using equation (2.3) for the intensity, integration of equation (3.3) over e and O yields for 
the absorption coefficient 

h r r ( e v 1_a ( r 



I r i^h) [if) °n(P)V-<x*0)d£d£l, (3.5) 

J Je min (9) \ m eC / \ n g/ 

where is the angle between dQ, and the direction of propagation of the 7-ray, and the domain for 
the f2 integration corresponds to the range of solid angles subtended by the disk. Equation (3.5) 
can be broken into two separate integrals by transforming from e to (3 using equation (3.2). The 
result is 

where the integration over j3 is expressed by the function 

(a) = a~ l f 2(3 (1 - /3 2 ) " 1 a 77 (/3) d/3 , (3.7) 
■/ 

plotted in Figure 1. 

Equation (3.6) can be used to calculate the absorption coefficient for a 7-ray propagating in 
an arbitrary direction at an arbitrary point in space. However, it is sufficient for our purposes to 
consider 7-rays propagating in a plane that includes the symmetry axis of the disk. Referring to 
the geometry indicated schematically in Figure 2, we assume that the 7-ray propagates in the (x, z) 
plane, where the z-axis is the symmetry axis of the disk, and we set the azimuthal angle if = 
along the x-axis. In general, we wish to relate dd to the differential area in the disk dA = R dR dip 
using 

cosr] RdRdif 

dn = , (3.8) 

where r\ is the angle between the X-ray trajectory and the normal to the disk surface, and I is the 
distance between the photon-photon interaction point and the X-ray emission point in the disk. 
Representing the location of the 7-ray using polar angle 5 and distance from the origin p, and 
assuming that the 7-ray trajectory makes an angle $ with the z-axis (see Fig. 2), our expression 
for the absorption coefficient becomes 

, n f _ Jo / E \ a T/ , f Ro f 2n , „,„, , / R \~ u cosridipRdR . . 
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where 

- 9 9 „ 9 „ _. . . . p cos(5 — <3?) — R cosG?sin<£ pcos5 

£ 2 = p 2 + i? 2 - 2 p R sin <5 cos </? , cos 6> = - '— , cos rj = '——^ — . 

(3.10) 

Note that a 77 oc E a , which verifies that softer X-ray spectra yield larger optical depths. Equa- 
tions (3.9) and (3.10) comprise our fundamental result for the 7 — 7 absorption coefficient for a 
7-ray of energy E propagating in the power-law X-ray field of an accretion disk. 

3.2. Absorption Coefficient for Propagation Along the Disk Axis 

One of the most interesting cases to consider is that of the attenuation of 7-rays propagating 
along the symmetry axis of the disk, since this scenario is of direct relevance for the interpretation 
of the 7-rays detected by EGRET from blazar sources. In this case equation (3.9) is simplified 
substantially due to the cylindrical symmetry, and the integration over tp is trivial. We can obtain 
the absorption coefficient for a 7-ray of energy E propagating outward along the symmetry axis 
at height z by setting p = z and 5 = $ = in equation (3.9) and transforming from R to 9 using 
tan# = R/z, yielding 

<l„ j f p \ a / z \ - u /-e max (z) 

a 1 JE,z) = ? J o- T m{a)[ — / (1 - cos 6) a+1 tan" w 9 sin 9 d6 , (3.11) 

c \2m e c 2 J \R g J 7e mm(z ) 

where 9 ma _^(z) [9 m i n (z)] is the maximum (minimum) value of the interaction angle 9 at height z. 
Neglecting gravitational curvature of the photon trajectories, we have 

cos9 - iz) = WTM ' cos ' min(z) = ' (3 ' 12) 

By writing equation (3.11), we have implicitly assumed that z 2> h, where h is the half-thickness 
of the disk. This condition is trivially satisfied for a thin disk, but the disk may also be considered 
to be thick, so long as we confine our attention to points sufficiently far above the surface. 

4. CALCULATION OF THE OPTICAL DEPTH 



We can use equations (3.9) and (3.10) to calculate the 7 — 7 optical depth r 77 for a 7-ray of 
energy E created at radius R and height z and propagating at angle $ with respect to the z-axis 
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by writing 

/>oo 

Tjj(E, R, z, <E>) = / a 71 (E, p, S, ®)d\, (4.1) 
■/ o 

where A is the distance traversed by the 7-ray since its creation, and p and 5 are given by 

9 „o ?,9 , x # + A sin $ z + A cos $ 
p = i? + z + A + 2 A (i? sin $ + 2: cos <3?) , sm<5 = , cos 5 = . 

P 

(4.2) 



In keeping with equation (3.9), we have assumed that the plane of propagation of the 7-ray includes 
the z-axis. 

4.1. Propagation Along the Disk Axis 

We can obtain some further simplification for 7-rays propagating along the z-axis, since then 
5 = <E> = R = 0, and p = z + A. The optical depth to pair production for a 7-ray of energy E 
created at height z above the center of the disk and propagating outward along the z-axis is then 
given by 

a 11 (E, z*) dz* , (4.3) 

where a 11 is given by equation (3.11). Unfortunately, one must resort to multidimensional numeri- 
cal integration to compute r 77 even in this highly symmetrical case. However, it is possible to obtain 
a closed-form expression for r 77 if z > Rq. Approximating the angular integral in equation (3.11) 
yields for the absorption coefficient 

(\ —2a — 4 / j-p \a 



where 

7T Iq er T fy(a) 



A = 



D \ 2q+4-w / D \ 2a+4-w 



RgJ \ Rg 



(4.5) 



(2a + 4 - u) c 

Substituting equation (4.4) into equation (4.3) then yields the approximate far-field solution for 
the optical depth 

, AR q ( z \- 2a ~ 3 ( E \ a 

valid if z > R . In § 5 we use this result to estimate the height of the 7-ray "photosphere," defined 
as the value of z such that r 77 = 1. We also calculate the 7 — 7 optical depth for the 7-ray blazar 
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3C 279 by focusing on three illustrative scenarios for the propagation of the 7-ray: (i) propagation 
along the disk axis; (ii) propagation along a line intersecting the disk axis; and (iii) propagation 
parallel to the disk axis. 



The observations of strong fluxes of X-rays and high-energy 7-rays from the blazar 3C 279 
in June 1991 allows us to place constraints on the location of the 7-ray emission site based on 7 — 7 
transparency requirements. Under the assumption that the contemporaneously observed X-rays 
were produced by a two-temperature accretion disk, we can use the expressions developed in §§ 3 
and 4 to calculate the optical depth traversed by the observed high-energy 7-rays. Since variability 
arguments (Makino et al. 1989) suggest that a significant fraction of the X-rays were probably 
produced in a relativistic jet, the results we will obtain for the optical depth r 77 are upper limits. 

The X-ray flux reported by Makino et al. (1993; Fig. 2) during the outburst had spectral 
index a = 0.68 and flux parameter Fq = 5.76 x 10~ 5 s -1 cm -2 (note that the flux given by their 
eq. [1] appears too large by a factor of 10). For this value of a, we obtain y = 1.6 (see eq. [2.1]) and 
*(a) = 0.245 (see Fig. 1), and for the redshift of 3C 279 (Z = 0.538), we obtain the luminosity 
distance (Lang 1980) D = Hq 1 c [Z + 0.5(1 - q ) Z 2 } = 5.65 x 10 27 cm assuming q = 0.5 and 
H = lOOkms -1 Mpc -1 . We can relate Fq to the disk-frame intensity parameter I using 



As an application, we consider first the simplest geometrical case, that of a 7-ray propagating 
away from the black hole along the axis of a two-temperature disk emitting with the power-law 
X-ray intensity given by equation (2.3). This scenario is useful for determining the minimum height 
above the disk at which the observed 7-radiation could have been produced in 3C 279 in June 1991. 
We can calculate the optical depth to pair production for a 7-ray of energy E created at height 



5. APPLICATION TO 3C 279 




1 -1 



(5.1) 



according to equations (2.8) and (2.9). 



5.1. Distance to the Gamma- Ray Emission Site 
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z above the center of the disk and propagating outward along the z-axis using equation (4.3). In 
order to calculate r 77 , we must assume values for the inner (R ms ) and outer (Ro) radii of the two- 
temperature, X-ray emitting region. In Figure 3 we set i? ms = 6 R g (accretion onto a Schwarzschild 
black hole) and Ro = 30 R g (a typical value for the outer radius; see SLE and EK), and plot 
log r 77 as a function of log (M/M Q ) and \og{z/R g ) for a 7-ray of energy E = 1 GeV using the X-ray 
data taken during the June 1991 EGRET flare of 3C 279. The optical depth exceeds unity in the 
cross-hatched area, and r 77 = 1 along the heavy solid line denoting the 7-ray photosphere, with 
height z = z-y-y. For reference, z = Ro along the upper dashed horizontal line, and z = R ms along 
the lower dashed line. The slope of the contours depends on how close to the disk the 7-rays are 
produced relative to R ms and Rq. For z < R ms , the contours become increasingly vertical because 
the absorption coefficient loses its dependence on z close to the disk. For z > Ro, the contours 
make a transition to the far- field power-law behavior given by equation (4.6), which reduces to 
r 77 ^ 2.4 x 10 16 (M/M Q )- 1 (z/ J R 3 )- 4 - 36 for the parameters used in Figure 3. 

In the case of 3C 279, we find that when M > 4 x 10 11 M Q , the intensity is so low that r 77 < 1 
regardless of the value of z. On the other hand, we see that observable 7-rays can be produced 
quite close to the black hole for reasonable values of M if z > R . For example, if M = 10 9 M , 
then r 77 < 1 for z > z 11 = 45 R g = 2 x 10~ 3 pc. This suggests that 7-rays produced at the top of 
the evacuated funnel at the center of a hot accretion disk (perhaps as a consequence of MHD waves 
generated in the surrounding hot-disk plasma; see Becker & Kafatos 1994), can escape without 
7 — 7 attenuation by the disk X-rays (although they may suffer attenuation by scattered X-rays on 
larger spatial scales; see Dermer Sz Schlickeiser 1994). 

The precise value of the outer radius Ro of the two-temperature region is uncertain, but 
it is likely to fall in the range 30 < Ro/R g ~ 100 (SLE; EK). In Figure 4 we therefore recompute 
the r 77 contours using the same parameters as in Figure 3, except Ro = 100 R g . Here again 
z = Ro and z = i? ms along the upper and lower dashed lines, respectively. The differences between 
Figures 3 and 4 are not dramatic, although the height of the 7-ray photosphere has increased 
from z 11 = 45 R g to z 11 = 70 R g for M = 10 9 M Q . This is due to the presence of X-ray photons 
emitted at large radii, for which the interaction angle 9 (and therefore the probability of pair 
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production) is enhanced. In this case, equation (4.6) for the approximate far-field solution reduces 
to r 77 = 3.6 x 10 17 (M/M o )- 1 (z/i? g )- 4 - 36 , in good agreement with Figure 4 for z > R . 

To investigate the effect of the angular dependence of the 7 — 7 absorption coefficient on 
the values obtained for r 77 , in Figure 5 we plot the height of the 7-ray photosphere, z^/R g , as a 
function of the black hole mass M/M Q for R ms = 6R g and uj = (uniformly-bright disk; dashed 
line) or uj = 3 (two-temperature disk; solid line). We compare results obtained using our two 
standard values for the outer radius, Rq = 30 R g {open squares) and Rq = 100 R g (open triangles). 
For uj = 3, increasing R from 30 R g to 100 R g causes only a modest increase in the height of the 
photosphere as mentioned above. However, going from uj = 3 to uj = while holding R constant 
results in a substantial increase in z 11 due to the increased amount of X-radiation impinging on the 
axially-propagating 7-rays from large radii, where the interation angle 9 has its maximum value. 
This effect is particularly strong for R = 100 R g . 

It is also interesting to consider the effect of varying the inner radius of the disk, -R ms , since 
this will also impact on the angular distribution of the X-rays. In Figure 6 we plot the photosphere 
height z 71 /R g as a function of M/M© for R ms /R g = 2, R /R g = (30,100), and uj = (0,3). 
These results may represent approximately situations involving disks around Kerr black holes, for 
which the radius of marginal stability can approach 1.2 R g (Shapiro & Teukolsky 1983). In this case 
equation (2.2) for the flux F(R) is no longer valid, although the flux still maintains the approximate 
power-law radial dependence F(R) oc R~ 3 (Eilek 1980). Obviously our model cannot be rigorously 
defended in this context since we have completely neglected relativistic effects, but, bearing this 
caveat in mind, it is interesting to note the modest decrease in the height of the photosphere going 
from Figure 5 to Figure 6 when uj = 3, due to the decrease in the flux-averaged value of the 
interaction angle 9. By contrast, the results depicted in Figures 5 and 6 agree closely when uj = 
(uniformly-bright disk), because then most of the emission comes from radii R 3> R ms , and the 
value of R ms is therefore irrelevant. 

In Table 1 we give values for the height of the 7-ray photosphere z 11 /R g , along with values 
for the associated light-crossing timescale (which may approximate the variability timescale) , t 11 = 
z-y-y/c, for E = 1 GeV, uj = 3, and various values of the black hole mass M/M@, the outer disk 
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radius Ro/R g , and the inner disk radius R ms /R g . Table 1 also includes values for the disk- frame 
intensity parameter I . Using our standard parameters M = 10 9 M Q , R = 30 R g , R ms = 6R g , 
and u> = 3, we find that i 77 ~ 2d, in agreement with the 7-ray variability timescale reported by 
Kniffen et al. (1993). 

5.2. Gamma-Ray Collimation 

The dependence of the 7 — 7 absorption probability on the angular distribution of the X- 
rays suggests that pair production may help to collimate the high-energy radiation escaping from 
blazars. We expect that for given values of the inner and outer disk radii R ms and Rq, the optical 
depth to pair production for 7-rays originating at a fixed location above the center of the disk will 
increase with increasing propagation angle <E>, mainly due to an increase in the flux-averaged value 
of the interaction angle 9, and to a lesser extent because the 7-rays remain close to the disk for a 
longer period of time. We shall continue to work with X-ray data obtained during the June 1991 
flare of 3C 279 in all of the examples. 

We can use equation (4.1) to calculate the optical depth r 77 for 7-rays of energy E created 
at a height z above the center of the disk and propagating at an angle <I> with respect to the z-axis. 
In order to explore the angular dependence of 7 — 7 attenuation, in Figures 7 and 8 we plot the 
transmission coefficient exp(— r 77 ) calculated using the 3C 279 X-ray data as a function of <I> for 
Rms = 6 R g , E = 1 GeV, M = 10 9 M Q , and various values of Ro and z. The transmission coefficient 
gives the fraction of the 7-rays that would escape to infinity in a given direction out of those emitted 
in that direction. We set Rq = 30 R g in Figure 7, Rq = 100 R g in Figure 8, and, to illustrate the 
importance of angular effects, we set to = 3 in Figures 7a and 8a, and u = in Figures 7b and 
8b. Strong "focusing" of the escaping 7-rays along the z-axis is apparent in all cases, and, for 
u = 3, the half-width of the angular distribution, A$, falls in the range 10° < A<I> < 20°. Similar 
results are obtained for A$ when co = (Figs. 7b and 86), but the overall level of attenuation is 
much greater in these cases due to the enhancement of the X-ray emission produced at large radii 
in the disk. In Table 1 we give values for A<E> evaluated at the photosphere, z = z 77 , for various 
model parameters. Our results for A<I> arc comparable to the observationally determined values 
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(e.g., Padovani & Urry 1992), and we therefore conclude that for acceptable disk parameters, 7 — 7 
attenuation alone is enough to produce the degree of 7-ray collimation implied by observations 
which associate extragalactic 7-ray emission with blazars. 

5.3. Shape of the Gamma-Ray Photosphere 

The high-energy 7-rays observed from many blazars are emitted nearly parallel to the axis 
of the jet, as is evidenced particularly by the common occurrence of apparent superluminal motion 
in the sources. It is therefore interesting to calculate, within the context of our model, the 7 — 7 
optical depth for a 7-ray created at radius R and height z, and subsequently propagating parallel to 
the 2-axis of the disk, which can be obtained by setting $ = in equation (4.1). It is particularly 
interesting to determine the shape of the 7-ray photosphere by varying the value of the 7-ray 
creation radius R. 

In Figure 9, we plot the photosphere height z 11 /R g calculated using the 3C 279 X-ray data 
as a function of R/R g for various values of M/M Q , assuming our standard parameters Rq = 30 R g , 
R ms = 6R g , and oj = 3. We find that the height of the photosphere increases dramatically as M 
decreases, and the shape of the photosphere becomes quite concave. Conversely, for large values 
of M, the photosphere becomes nearly flat. We interpret this behavior as follows. For large values 
of M, the disk-frame intensity parameter I is small, and therefore the photosphere is located 
relatively close to the disk, with z 77 /i?o < 1. In this case, most of the 7 — 7 absorption takes place 
so close to the disk that the flux-averaged value of the interaction angle 9 is nearly 90° at all points 
in the photosphere, and therefore the height of the photosphere varies little across the disk. On the 
other hand, as M decreases, Iq becomes large, and consequently the photosphere lies further from 
the disk, with z 11 /Rq > 1. In this case, the flux-averaged value of 9 increases significantly going 
from the center of the photosphere to the edge, and therefore the height of the photosphere must 
increase rapidly in order to keep r 77 = 1, as seen in Figure 9. 

6. DISCUSSION 

Although the physics of particle creation and acceleration far from the turbulent plasma 
accreting onto a black hole is poorly understood, there are a number of relatively well-understood 
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plasma and accretion instabilities that can give rise to efficient particle acceleration within the 
disk itself, or within the evacuated funnel region at the center of the disk. For example, resonant 
turbulent acceleration via interactions between particles and (linear or nonlinear) Alfven waves can 
in principle lead to the production of electrons with very high Lorentz factors, and to the associated 
emission of high-energy inverse-Compton 7-rays (e.g., Becker & Kafatos 1994; Becker, Kafatos, k, 
Maisack 1994). However, in order for any of these mechanisms to create observable 7-rays, the 
optical depth to pair production on the background UV and X-ray photons must be small enough 
for the 7-rays to escape from the source. 

In this paper, we have attempted to bridge the gap between theory and observation by 
developing useful constraints based on the simple hypothesis that the power-law X-ray spectrum 
originates in the inner region of a two-temperature accretion disk. Such disks commonly provide 
good fits to the hard X-ray spectra of AGNs. In § 3 we developed a general expression for the 
7 — 7 absorption coefficient a 77 for 7-rays propagating above a two-temperature disk in a plane 
that includes the symmetry axis, and in § 4 we derived an approximate far-field solution for the 
7 — 7 optical depth r 77 for 7-rays created above the center of the disk and propagating outward 
along the disk axis, taking full account of the angular dependence of the X-ray distribution. In § 5 
we applied our results to 3C 279 using the X-ray data reported by Makino et al. (1993) during 
the June 1991 7-ray flare, and demonstrated that the observed 7-rays could have been produced 
within < 45 GM/c 2 of a black hole of mass M > 10 9 M Q . This is close enough to the central source 
to establish a direct connection between the production of the 7-rays and the accretion of material 
onto the black hole. In particular, our results suggest that the observed 7-rays may have been 
produced near the top of the evacuated funnel at the center of a hot, two-temperature accretion 
disk. In this scenario, the relativistic electrons that Compton-upscatter UV photons to produce the 
7-rays are themselves produced as a consequence of photon-photon collisions in the funnel. The 
absence of protons in the funnel due to the centrifugal barrier allows the presence of Alfven waves 
with wavelengths short enough to effectively accelerate the electrons (Becker, Kafatos, Sz Maisack 
1994). 

In order to explore the observational implications of the angular dependence of 7 — 7 pair 
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production, in § 5 we also computed r 77 for 7-rays born along the disk axis but propagating at 
an angle with respect to it. The results for the half-angle (A<I> < 15°) indicated in Figures 7 
and 8 and in Table 1 suggest that pair production in the X-ray field of the accretion disk may 
be strong enough to explain the observed degree of 7-ray alignment in blazars. Furthermore, the 
variability (light-crossing) timescale i 77 = z 11 /c we derive for black hole mass M ~ 10 9 M Q is in 
good agreement with the 2-day 7-ray variability of 3C 279 reported by Kniffen et al. (1993). 

The results presented in §§ 4 and 5 suggest that 7 — 7 pair production can help to "focus" an 
intrinsically isotropic source of 7-rays so that the high-energy radiation escapes preferentially along 
the disk axis. Whether or not the 7-rays actually are produced isotropically is a separate question. 
If they are, then several additional considerations come into play, regarding (i) the self-opacity of 
the 7-rays; (ii) the fate of the 7-ray energy emitted off-axis; and (iii) the required total luminosity. 

6.1. 7 — 7 Self-Opacity 

Whether or not the high-energy 7-rays are produced by an intrinsically isotropic source, their 
self-interaction optical depth r se if cannot exceed unity if they are to be observed. The most severe 
self-interaction constraint is obtained in the case of isotropic emission. Assuming a monochromatic, 
isotropic distribution of 7-rays with energy Eq > m e c 2 and intensity I e = B S(E — Eq), integration 
of equation (3.3) yields for the self-interaction absorption coefficient 
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where cr 77 (/3) is given by equation (3.1) and 
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If the 7-ray source is an isotropically emitting sphere of radius i?* , then the intensity parameter B 
is related to the 7-ray luminosity L 7 via 



(2vri?*) 



B = ToTT^I > ( 6 - 3 ) 



and we obtain for the self-interaction 7 — 7 optical depth 
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An isotropic 7-ray source can be ruled out if r se if exceeds unity for reasonable values of R*. The first 
factor on the right-hand side of equation (6.4) is essentially the standard 7-ray compactness, and 
the quantity in braces becomes ~ 10~ 6 for E = 1 GeV. Taking L 7 = 10 48 ergss _1 for 3C 279 and 
setting T S eif < 1 yields the transparency requirement i?* > 10 10 cm, and therefore self-interaction 
of the GeV radiation in 3C 279 is not likely to be important. A similar calculation by McNaron- 
Brown et al. (1995) indicates that the self-interaction of radiation in the energy range 0.05 — 10 MeV 
(observed by OSSE) is also negligible. 

6.2. Fate of the Off-Axis Radiation 

If the 7-rays are produced isotropically, then only those emitted close to the disk axis will 
escape from the X-ray field of the accretion disk. The radiation emitted further from the axis will 
be degraded in a pair-Compton cascade (Protheroe, Mastichiadis, & Dermer 1992), with most of 
the luminosity residing ultimately in a broad MeV feature. In this view, objects such as 3C 279 
(while in 7-ray flare mode) will appear as bright sources of MeV emission when viewed along lines 
of sight > 20° from the disk axis. Centaurus A may be an example of such an object (Skibo, 
Dermer, & Kinzer 1995). 

The luminosity required to power isotropic 7-ray emission, L 7 = 10 48 ergss _1 , is orders 
of magnitude higher than that required if the emission is beamed. This is perhaps the simplest 
argument against isotropic emission, but it is especially compelling when coupled with the strong 
circumstantial evidence in favor of relativistic bulk motion. However, there is apparently no strong a 
priori reason to rule out isotropic emission, since attenuation due to collisions with X-rays emitted 
by the disk is sufficient to produce the observed degree of alignment. In any event, our results 
regarding 7 — 7 collimation do not rely on the existence of an isotropic 7-ray source. The mechanism 
we consider here can work in concert with other sources of collimation, such as relativistic bulk 
motion or geometrical focusing. 

6.3. Attenuation by Scattered X-Rays 

The backscattering of X-rays produced in the accretion disk by electrons located above the 
disk tends to isotropize the X-ray distribution, and this can have a strong effect on the results 
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obtained for the 7 — 7 absorption coefficient. We believe that our neglect of a scattered component 
to the X-ray distribution is reasonable due to the proximity of the 7-ray production site to the black 
hole in our model. It has been shown (Dermer k, Schlickeiser 1994) that accretion disk photons will 
dominate over scattered X-rays within a distance of ~ 0.01 — 0.1 pc from the central source, which 
is much larger than the typical height of the 7-ray photosphere, z 77 , obtained here (see Table 1). 
However, the problem of 7 — 7 attenuation by scattered X-rays on much larger spatial scales is 
a serious issue that must be addressed using detailed calculations that determine self-consistently 
the angle and energy distributions of both the scattered radiation and the electron-positron pairs 
at large distances from the black hole. Such a calculation is beyond the scope of this paper, but it 
suggests a promising direction for future work. 

Our results indicate that if the emission height z is held fixed, then the optical depth increases 
with increasing 7-ray energy E as r 77 oc E a , and therefore the height of the 7-ray photosphere z 77 
increases with increasing E. If we arbitrarily set the emission height for all of the 7-rays equal to 
the value of z 11 obtained for E = 1 GeV, then the observed spectrum will be attenuated by a factor 
~ 7 going from E = 1 GeV to E = 5 GeV. It is unclear whether such a steepening actually exists 
in the spectrum of 3C 279 (Hartman et al. 1992). Finally, for energies E > 5 GeV, the interaction 
with UV radiation produced in the cool outer region of the disk would also become important, 
further raising the height of the photosphere. 

To summarize our main points, application of our general expressions describing 7 — 7 
attenuation above X-ray emitting accretion disks leads us to conclude that in the case of 3C 279, 
(i) 7-ray transparency (to disk-generated X-rays) is achieved if the 7-rays are produced at height 
z > Ro; (ii) 7 — 7 attenuation is strong enough to explain the observed association of high-energy 
7-ray emission from AGNs exclusively with blazars, and (iii) the observed 7-rays may be produced 
in active plasma located above the funnel at the center of a two-temperature accretion disk. 

We gratefully acknowledge the remarks provided by the anonymous referee, whose insightful 
comments contributed significantly to the paper. This work was partially supported by a NASA 
Cycle 3 CGRO Guest Investigator grant. 
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FIGURE CAPTIONS 

Fig. 1. - The cross section integral function *(a) (eq. [3.7]). For the 3C 279 data, a = 0.68 (see 
Makino et al. 1993), and ^ = 0.245. 

Fig. 2. - Schematic representation of the 7-ray propagation geometry. The 7-ray is located at 
distance p from the origin and polar angle 5 and propagates in a plane that includes the z-axis. 
Attenuation of the 7-rays occurs as a result of collisions with X-rays generated at all points 
on the disk. The interaction angle between the photons is 9, and the angle between the 7-ray 
trajectory and the z-axis is <E>. 

Fig. 3. - Contour plot of the log of the 7 — 7 optical depth r 77 (eq. [4.3]) for 7-rays with energy 
E = 1 GeV created at height z above the center of a two-temperature disk (lo = 3) with outer 
radius R = 30 R g and inner radius R ms = QR g . The upper and lower dashed horizontal 
lines denote z = Ro and z = R ms , respectively, and the 7-rays propagate outward along the 
z-axis. The curves were computed using the X-ray data for 3C 279 taken during the June 1991 
EGRET flare by Makino et al. (1993). 

Fig. 4. - Same as Fig. 3, except R = 100 R g . 

Fig. 5. - Height of the 7-ray photosphere, z 11 /R g , plotted as a function of the black hole mass 
M/Mq for i? ms = QR g and lo = (dashed line), lo = 3 (solid line), Ro = 30 R g (open squares), 
Ro = 100 R g (open triangles). The curves were computed using the June 1991 3C 279 X-ray 
data. 

Fig. 6. - Same as Fig. 5, except R ms = 2R g . 

Fig. 7. - Transmission coefficient exp(— r 77 ) plotted as a function of the propagation angle using 
the June 1991 3C 279 X-ray data for M = 10 9 M , R = 30 R g , R ms = 6R g and (a) uj = 3, 
(b) lo = 0. The 7-rays are created at the height z/R g indicated for each curve. 

Fig. 8. - Same as Fig. 7, except Ro = 100 R g . 

Fig. 9. - Height of the 7-ray photosphere z 11 /R g plotted as a function of the creation radius R/R g 
for 7-rays propagating parallel to the z-axis using the June 1991 3C 279 X-ray data, assuming 
Ro = 30 R g , R ms = QRg, and lo = 3. The value of the black hole mass M/M Q is indicated for 
each curve. 



